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NCNUNIFORMITIES IN SHOCK- TUBE FLOW DUE TO 
UNSTEADY-BOUNDARY-LAYER ACTION 
By Harold Mirels and W. H. Braun 


SUMMARY 

The boundary layer along the walls of a shock tube induces pressure 
and velocity gradients within the core of potential flow. These nonuni- 
formities are evaluated herein for shock tubes in which the boundary 
layer is thin relative to the tube diameter and is either wholly lamina r 
or wholly turbulent . The hot gas region between the shock wave and the 
contact surface is considered. Both the axial distributions at any in- 
stant and the temporal distributions at any axial position within this 
region are found. Numerical computations are presented for an air-air 
shock tube . 


INTRODUCTION 

The analysis of reference 1 obtains the flow pert urbat ions, due to 
unsteady-boundary-layer action, immediately behind the shock wave in a 
shock tube. In this report the method will be extended to find the flow 
perturbations in the entire region between the shock wave and the con- 
tact surface. These perturbations are of interest when experimental 
shock- tube data are analyzed. As in the previous work, the assumptions 
are made that the boundary-layer action is equivalent to a one- 
dimensional distribution of mass sources and that the expansion wave has 
zero thickness . Numerical results are obtained for an air-air shock 
tube for which the initial temperature throughout the tube is 520° R, 
the wall boundary layer is either wholly laminar or wholly turbulent, 
and the shock- tube wall remains at a temperature of 520° R. The spatial 
and temporal variations in the hot gas region are presented graphically. 

A theoretical study of the perturbations in the hot gas region be- 
tween the shock and the contact surface is reported in reference 2 . The 
limitations of the method of reference 2 are discussed in reference 1. 

A comparison of some of the results of reference 2 with those of the 
present report is made in the body of this report . It might also be 
noted that the method of reference 2 is also applied in reference 3 to 
estimate the perturbations in the expansion fan of a shock tube. 
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ANALYSIS 
Perturb at ions 

In the following analysis, familiarity with reference 1 is assumed. 
The region between the shock wave and the contact surface is considered. 
The perturbation at an arbitrary point in this region (point d(x,t) in 
fig. l) is found by summing the pressure impulses which arrive at that 
point from both the left and the right. This summation requires an in- 
tegration along all those characteristics influencing point d(x,t). The 
characteristic lines considered herein are indicated in figure 1. Ex- 
pressions for these characteristic lines sire given by the following equa- 
tions in terms of £,t. The equations for u g x - |, included therein, 

will be useful in later developments in this report. (Symbols are defined 
in appendix A . ) 
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Following the method of reference 1, the contribution to the pres- 
sure perturbation at point d in region 2 (hot gas) from waves arriving 
from the left is 
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The integrand is, in each case, the vertical velocity at the edge of the 
boundary layer. Coefficients C and D are, respectively, the reflec- 
tion coefficient and the transmission coefficient at the interface between 
regions 2 and 3. They are derived in reference 1, appendix C, and equal: 


'\ 


C = 


D = 


P 32 a 23 


- 1 


T 32 a 23 + 1 




( 3 ) 


r 32 a 23 + 1 


y 


4 


MCA TN 4021 


The waves arriving at point d from the right contribute a pressure 
perturbation 
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The reflection coefficient E at the shock is defined in reference 1, 
appendix C , as 
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The contributions to Ap 2 d and Ap 2 d from waves reflected at points 
c and c are shown to be negligible in reference 1. 

The limits of the integrals in equations (2) and (4) are the end 
points of the characteristic line segments (eqs. (l)). They may be 
evaluated in terms of the coordinates x,t at point d as follows: 
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The vertical velocity at the edge of the boundary layer is ex- 
pressed in reference 1, appendixes D and E, as 



The two coefficients L 2 and L 3 are evaluated in reference 1 for both 
wholly laminar and wholly turbulent boundary layers. The indices n 2 
and n 3 are l/2 for laminar boundary layers and l/5 for turbulent 
boundary layers. 
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The use of equations (6) and (7) leads to the following values of 
the integrals in equation (2): 
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Under tbe assuinptions that "the boundary layer is wholly laminar or 
wholly turbulent (n 2 = n 3 = n), the two pressure disturbances become 
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The complete pressure and velocity perturbat ions at point d are 
obtained from equation (9) according to 
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Other perturbation quantities, such as temperature and density, can be 
found by the methods indicated in appendix B. 

The axial distributions of the perturbations Ap 0 , and Au„ , at 

^ y U. Cj y Q. 

any instant of time are presented in figures 2 and 3 for several nominal 
shock numbers . An air-air shock tube with = T^ = 520° R and wholly 

laminar or wholly turbulent boundary layers is considered. The wall tem- 
perature is assumed to remain constant at 520° R at all times. The ratio 
of specific heats x and the Prandtl number were taken as 1.4 and 0.70, 
respectively. In each case the curve is normalized to the value Apg 0 

(or Aug Q ), which is the value of the perturbation immediately behind 

the shock wave at the instant of time under consideration. In using the 
plots of Ap_ -|/Ap P , for example, the numerical value of Ap_ may 

be obtained from equations (9) and (10a) evaluated at x = u s t. A non- 
dimensional plot of Ap 2 q is given in reference 1 for an air-air shock 

tube and is reproduced here for convenience in figure 4. If an experi- 
mental observation of shock attenuation is made, Ap_ can be determined 

Cl y O 

directly from the experimental data (see, e.g., appendix B). Then fig- 
ure 2 should yield accurate estimates of the axial variation of pressure 
between the shock and the contact surface at the fixed instant of time. 

If equations (9a) and (9b) are multiplied by (u s t/x)^“ n , they are 
brought into a form which gives the variation, with time, of the pertur- 
bations at a fixed axial station x (appendix B). This form is of in- 
terest when aerodynamic measurements are made at a fixed station in the 
shock tube. Numerical results are presented in figures 5 and 6 (for the 
air-air shock tube described in connection with figs. 2 and 3). The 
abscissa is the dimensionless time measured from the passage of the shock 
past the fixed station under consideration. The curves terminate at that 
time at which the contact surface passes the station. In figures 5 and 
6, Apg 0 is the pressure perturbation directly behind the shock at the 

instant it passed the fixed station being considered. The numerical 
value of Apg 0 can a S a ^- n be determined by the methods described in 

appendix B. 

Typical results of the present analysis are compared with those of 
reference 2 in figure 7 . The temporal variation of pressure at a fixed 
station in an air-air shock tube is shown therein. The temperatures are 
^ = T 4 = 520° R; the axial position is 8 feet and the hydraulic diameter 

d = l/7. The solid curves are based on equations (9). The dashed 
curves are from figure 8 of reference 2. It is seen that, for the range 
of Mach numbers considered, the method of reference 2 predicts a greater 
increment of the perturbation pressure in the hot gas than does the 
present method. For example, consider conditions at the contact surface 
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for p 21 = 3.4. The ratio Ap 2 ^/Apg o e Q. ua l s -0-2 by the method of 

reference 2 and equals 0.2 by the present method. This represents a 
discrepancy of about 40 percent in the perturbation quantity Apg (j/Apg^ 

The limitations of the method of reference 2 are discussed in reference 


Weak: Shocks 


If the same gas is used in regions 1 and 4 = 14 ) at the same 

initial temperature (T^ = ) , and if the boundary layers are wholly 

turbulent or wholly laminar (ng = ) , then, for weak shocks (M g = 1, 

Mg « l) the pressure and velocity perturbations may be found in terms 
of the single parameter Mg. 


By using the relation 
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which may be obtained from the ideal shock- tube relations (as given, 
e.g., in appendix G of ref. l), the transmission and reflection coeffi- 
cients are found to be 
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By neglecting higher-order terms in Mg, equations (9a) and (9b) become, 
respectively, 



valid for Mg < x/u g t < 1. (The point x/u g t = Mg corresponds to the 

contact surface.) Equations (ll) follow from equations (9) by neglect- 
ing the term containing C in equation (9a) and the terms containing E 
in equation (9b). The pressure and velocity perturbations in region 2 
(eqs. (lOa) and (b), respectively) are now found to be 



The limiting forms (12a) and (l2b) are shown in figures 8 and 9, 
respectively, for three values of Mg. The accuracy of the limiting 

forms is indicated by comparison with the dashed curves, which were ob- 
tained from equations (9) and (10) for Mg = 0.2 (M g = 1.1328). The 

corresponding time plots are shown in figures 10 and 11. They are ob- 

1 — n 

tained as the products of equations (12) and (u g t/x) 
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DISCUSSION OF RESULTS 


The numerical results of the present analysis are presented in 


graphic form in figures 2 to 11. The normalized ratios 

and Au„ ,/Au„ are used therein. 

2 ^ ci. 2, o 


2,d / "" ± '2,o 
It should he kept in mind that 


are negative in all cases so that when these normal- 


AP 2 ,o and Au 2 ^ 0 

ized ratios are positive or negative, the numerator is negative or posi- 


c pc 

tive, respectively. The ratios Ap 2 &fa?2,o 
and Au 2 d /Au 2 0 are less when variations 

with x at a fixed t are considered than 
when variations with t at a fixed x (ex- 
cept at x/u g t = l) are considered, since 

Ap 0 and Au„ are larger for the former 
2,o 2,o 

than for the latter (see sketch) . The re- 
sults for wholly laminar boundary layers ex- 
hibit the same trends as those for wholly 
turbulent boundary layers, although the individual perturbations are less 
for the laminar case (ref. l) . Hence, only the trends of the turbulent 
boundary -layer results will be discussed herein. 



For weak shocks the ratio Ap 2 ^/Ap 2 Q varies from 1 at the shock 

to nearly zero at the contact surface (figs. 8(b) and 10(b)). Thus, 
the pressure near the contact surface remains nearly at the ideal (i.e., 
no-attenuation) value. The ratio Au^ ^/Au^ 0 varies from 1 at the 

shock to a substantial negative value at the contact surface (figs. 9(b) 
and 11(b)). Thus the velocity perturbation at the contact surface is 
positive and the contact surface "speeds up" because of the wall boundary 
layer. This is observed experimentally in reference 4. 

With increases in shock strength, the value of Ap 2 ^/Ap 2 ^ Q at the 

contact surface also increases (figs. 2(b) and 5(b)) so that a consider- 
able pressure decrement exists at the contact surface . For strong shocks 
(Mg = 0(6)) the pressure perturbation at the contact surface is of the 
same order of magnitude as that behind the shock so that the pressure is 
fairly uniform in region 2. This result is in qualitative agreement with 
the unpublished experimental measurements of J. J. Jones of the Langley 
laboratory, in which a relatively constant pressure region directly be- 
hind the shock is noted. With increase of Mg, the value of Au 2 ^ d /Au 2 ^ 0 
at the contact surface ultimately increases to a value of about 0 for 
Mg * 6 (figs. 3(b) and 6(b)). Thus for Mg * 6 the contact surface 
moves at approximately the ideal value. However, since the shock speed 
is lower than ideal, due to attenuation, the speed of the contact surface 
relative to the shock is greater than that predicted by ideal shock-tube 
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theory. This effect tends to reduce the available test time when ex- 
periments axe conducted at a fixed station in a shock tube. Moreover, 
this results in a substantial variation of velocity in region 2 despite 
the fact that the pressure variation is relatively small. 

The previously noted results can be explained on the following 
physical grounds. The boundary layer in region 3 generates longitudinal 
compression waves while that in region 2 generates expansion waves (ref. 
l) . The net effect of these waves is to attenuate the shock as described 
in reference 1. The shock attenuation defines the perturbations at one 
limit (x/u s t » l) of region 2. The variation of fluid properties between 

the shock and the contact surface can then be estimated (approx.) by con- 
sideration of conditions near the contact surface (x/ugt ~ Ug/ug). 


For weak shocks the perturbations near the_contact surface are 
mainly due to the characteristic lines ab and dd (see sketch) . The in- 
tegration along line ab contributes, at 
point d, a compression wave moving in the 
+x-direction. The integration along line 
dd contributes an expansion wave moving in 
the -x-direction. For an air-air shock 
tube with small initial pressure ratio, 
these pressure waves are approximately of 
equal and opposite strength so that the 
pressure perturbations cancel each other 
(eq. (10a)) while the corresponding veloc- 
ity perturbations are additive (eq. (lOb)). 
Thus for weak shocks the pressure pertur- 
bation is essentially zero at the contact surface while the velocity 
perturbation has a substantial positive value. Other perturbation quan- 
tities can be readily found by using appendix B. 



Now, consider the case of an air -air shock tube with a very high 
initial pressure ratio so that the shock strength is large. It can be 
shown that region 3 (i.e., characteristic line ab) makes only a s m all 


Region 3 


^ Region 2 



contribution to the pertur- 
bations at points near the 
contact surface. This is 
primarily due to the fact 
that the transmission coeffi- 
cient D at point b (i.e., 

D = ^>2,b/^3,b iS the ra ‘ 
tio of the wave transmitted 
into region 2 to the inci- 
dent wave arriving from 
region 3) approaches zero 
for the strong- shock case. 
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Thus the major contribution to the perturbations is due to characteristic 
lines dd and cb. The integrations of Vg along lines dd and cb are very 

nearly equal when point d is near the contact surface. However, the con- 
tribution of line cb to point d is found by multiplying the integral 

along line cb by the reflection coefficient C (i.e., C = Apt , /An~ , is 

’ 2,b' 2,D 

the ratio of the reflected to the incident wave, in region 2, at point 
b). But C approaches 1 for strong shocks so that the pressure pertur- 
bation_at point d, due to line cb, is very nearly equal to that due to 
line dd. However, since the contribution of line cb is a wave moving 
in the +x-direction (due to reflection at point b) while line dd contri- 
butes an approximately equal pressure wave moving in the -x-direction, 
the pressures are additive, but the resulting velocity perturbation is 
nearly zero. As a result the pressure decrement at the contact surface 
is of the same order as that behind the shock, but the velocity of the 
contact surface remains at essentially the ideal value. (The latter 
result could also be deduced by noting that the acoustic impedance ratio 
P 2 a j/p 2 a 2 = ^3 a 2A"2 a 3 as ver y s0 that the contact surface tends to 

be a constant-velocity interface and at the pressure imposed by region 2.) 


The perturbations at a fixed station x for the case of a strong 
shock is discussed further in appendix B. It is noted therein that, 
during the time interval between passage of the shock and passage of the 
contact surface, Apg ^/Apg 0 remains at about 1, Aug ^/Aug ^ Q varies 


from 1 to 0, 


Ap 


2 ^ 8 . 


P2 



varies from 0 to 


1_ ^~2 ,d p 2 
*l' T 2 Ap 2,o 


varies froml 


T-, - 1 ^2 d P2 Ti - ! 

to — = , and ■ — varies from 0 to — . In interpret- 

r l ^ Ap 2,o ^1 

ing these results, it should be noted that Ap and Au_ are neg- 

Lx yO Lx y O 

ative. Thus, there is only a slight increase in AMg^, and Apg^ re- 
mains essentially constant, during the course of an experiment employing 
an air-air shock tube at large pressure ratios. 


CONCLUDING REMARKS 

The calculations which have been made herein are of use in several 
kinds of shock- tube studies. Probably, the most important application 
is the prediction of the variation of free- stream conditions at a sta- 
tionary model (or instrument) utilizing the hot gas region of a shock 
tube. If the shock tube is to be used for the study of boundary-layer 
stability or dissociation effects on the development of the boundary 
layer behind the shock, then the nonuniformities derived herein are also 
of interest . 
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The limitations of the present method are the same as those dis- 
cussed in reference 1. The use of a zero-thickness expansion wave is 
valid for small pressure ratios but becomes in error for large pressure 
ratios. However, for very large pressure ratios the contributions of 
region 3 (to conditions in region 2) are negligible due to the acoustic- 
impedance mismatch at the contact surface. Hence, it is only for inter- 
mediate shock strengths that the use of a zero-thickness expansion wave 
is likely to introduce errors . The assumption that the boundary layer 
generates only longitudinal waves introduces a steep gradient of the 
perturbation quantities directly behind the shock (e.g., fig. 2). (This 
is due both to the singular nature of the boundary layer and to the 
assumption of one -dimensional flow within the potential-flow core.) The 

gradient is steeper for the laminar case (when Vg has a jj - order singu- 
larity directly behind the shock) than for the turbulent case (where Vg 
has a —order singularity). If the three-dimensional nature of the flow 

is taken into account, the steep gradient behind the shock would be ex- 
pected to be somewhat modified. In experimental studies of the relax- 
ation effects associated with strong shocks, it may be necessary to dis- 
tinguish between the perturbations due to relaxation effects and those 
due to the wave system introduced by the wall boundary layer. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 12, 1957 
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APPENDIX A 
SIMBOLS 

speed of sound 

reflection coefficient at interface, eq. (3) 

specific heat at constant volume 

transmission coefficient at interface, eq. (3) 

hydraulic diameter, ft 

reflection coefficient at shock, eq. (5) 

eq. (7); ref. 1, appendixes D and E 

Mach number of flow relative to wall 

Mach number of shock wave relative to wall 

eq. (7); ref. 1, appendixes D and E 

pressure 

temperature, °R 

time, sec 

velocity of shock wave relative to wall 
velocity in region 2 relative to wall 
vertical velocity at edge of boundary layer 
longitudinal distance from diaphragm 
ratio of specific heats 

perturbation from ideal (no attenuation) shock- tube flow 
kinematic viscosity 

integration variable representing x 
mass density 
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T 

• • • 
a,b, • • • 

Subscripts: 

points in x,t diagram (fig. l); subscript 2,d refers to 
arbitrary point in region 2 

condition at shock wave; subscript 2,0 refers to conditions in 
region 2 directly behind shock 

wall 

regions of shock tube (fig. l) 

Superscripts: 

+, - associated with waves moving in + or - x-direct ion, 

respectively 

Special notation: 

a „ = a,/a^; L™ = L,/l» * * • * two successive integer subscripts, not 
43 4' 3 32 3' 2 

separated by a comma, represent a ratio 


a ,b, • • • 
a,b, • • • 

o 

w 

1 , 2 , 3, 4 


integration variable representing t 

points on characteristic lines (fig. l); point d is an arbitrary 
point in region 2 
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APPENDIX B 


PERTURBATIONS IN REGION 2 

Useful expressions relating the perturbation quantities in region 
2 are derived herein. The perturbations represent the departure from 
the ideal conditions that would exist in region 2 if the wall boundary 
layer were not present. A perfect gas is assumed so that y-g = T]_* 


Perturbations Directly Behind Shock 


Under ideal conditions, regions 1 and 2 are related by the normal- 
shock relations 


p 2 2r ± 4 - (r x - i) 


Pi 

Ti * 1 

^2 

2 4 - 1 

a l 

" r 1 + 1 Mg 

£2 

(rp + i)m| 


(bi) 


p i ( r - i)m| + 2 


j 


where M is the shock Mach number. If the shock is perturbed, con- 
s 

ditions directly behind the shock (designated by subscript 2,o) can be 
expressed as 
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where As 


4 Ti m2 


^s 


^2,o = _ 

p 2 2r x M| - (r x - i) **8 

^s - frl - Ap 2,( 
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4ri^ 
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4 ^i M s 




1 P2 
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^2,0 _ ^l M s - Cri - 1) Ap 2>Q 

p 2 rM 2 [(ri - 1)M? + 2 ] p 2 

S -L ^ 

^-2 = — ^-2 ^-2 (Equation of state) 


^2,0 _ ^ p 2,o ^2,0 

c v p 2 ~ T p 2 

frl - 3-)(nj - if 4?2 ,o 

m l[{r ± - i)m| + 2] P2 

is the entropy perturbation. 
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Entropy Variation in Region 2 

The entropy is constant for a given particle in region 2, the value 
depending on the strength of the shock at the instant the particle passed 

through it. A particle at x,t in 
x region 2 passed through the shock 
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The entropy of a particle at x,t is then given by 

^2, a frl - 1 )( m | - J -) 2 ApgA 

° v ■ ^ICr, - + *1 V** 4, tA CB4) 

For weal; shocks, equation (B4) is of the order of (M - l) 3 and the en- 
tropy variations are negligible. The entropy perturbations do not in- 
teract with the longitudinal pressure waves, to the order of the present 
analysis, so that each can be considered separately. 


Parametric Form, of Perturbation Solutions 


For wholly laminar boundary layers (n = l/2) or wholly turbulent 
boundary layers (n = l/5) throughout the shock tube, the local pertur- 
bation of any quantity in region 2 (designated by A( ) 2 can be ex- 
pressed in the fonns 



\ 




(B5) 


where 




The form of equation (B5) can be de- 


duced from equations (9) . The function f| 
the variation with x at a fixed time t 


fe) 


can be interpreted as 
while the function g 


■fe) 


can 


be inteipreted as the variation with t at a fixed x. Both forms are 
used in figures 2 to 11. When dealing with perturbations normalized to 
conditions directly behind the shock (i.e., A( ) 2 d /A( ) 2 Q ), multipli- 
cation by (u s t/x) 1-n converts a solution for fixed t to a solution 
for fixed x. 


u s t u s 

Perturbations at Fixed x for Times 1 < ^ — 

X Ug 

Let subscript 2,o be considered to represent the perturbation 
directly behind the shock at the instant it passes a fixed station x 
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and let subscript 2,d represent the 
perturbation at the station at a sub- 
sequent time. Conditions directly 
behind the shock can be found by (a) 
experimental determination of ZMg 

or Apg^ 0 /P2 equation (B2) or 

(b) analytical determination of 
Ap 2 o/ p 2 f ’ rom equations (9) (or fig. 

4 for an air-air shock tube) and (B2) 
The subsequent perturbations can then 
be expressed as indicated in the fol- 
lowing section. 


The entropy perturbation may be written in the form 


^ 2,1 p 2 

c v Ap 2,o 


( Tl - 1)(M^ - i) 2 

M g[(Ti - l)M| + 2] 



(B6) 


which is obtained from equations (B4), (B5), and (B3) by noting 



Other perturbations of interest are 


^2,d p 2 = 1 f ^2,d ^2 ,d p 2 \ 

p 2 Ap 2,o = r \^2,o " C v Ap 2,oJ 


(B7a) 


A ^2,d p 2 ^1 ~ 1 Ap 2 ,d | ^2^ p 2 

T 2 Ap 2,o " Ti AP 2 ,o + Ti c v Ap^ 0 


(B7b) 


^2,d p 2 
•^2 ^ 2,0 


Au 2 ,o p 2 \ Au 2,d 
v «2 Ap 2,oy ^ 2,0 


1 ( m 2,d p 2 \ 

2 V T 2 Ap 2 ,oJ 
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^2,(1 p 2 _ ^2,d p 2 + 2 f ^>0 p 2 \ ^ u 2,d 

«2 Ap 2,o p 2 Ap 2,o + \ u 2 Ap 2, 0 y ^ 2,0 


(B7d) 


where Mg = Ug / a 2 and q 2 = P 2 u|/2. 1116 ratios Ap 2 d/^ p 2 o and - 
Au 2,d/ Au 2,o are found from equations (9) and (10) or from figures such 


Au< 


as 5 and 6. The ratio 


2^0 


p 2 


"2 ^ 2,0 


can be found from equation (B2) 


It is of interest to investigate the nature 

for very strong shocks . Assuming » 1 and 
tions (B2) show 


of the perturbations 


(ri - i)^ 

P » 1, equa- 
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(B8) 


which defines the perturbations directly behind the shock. The subse- 
quent time variation of the perturbation can then be estimated by con- 

u s t u s 

sidering conditions at the contact surface, = — . For an air-air 

x u 2 

shock tube and Mg * 0(6), figures 5 and 6 indicate that at the contact 
surface 


Au 

Au 


2,d 

2,0 


^2,d 

2,0 


0 

> 

1 

J 


(B9) 


Equations (B9) probably are applicable for other driver gases provided 
Mg is sufficiently large (in order to make p 2 a 2 lar se) so that 

region 3 has only a small effect on region 2. (See section titled 
Discussion of Results.) By use of equations (B9), the other perturba- 
tions at the contact surface are 


^2 ,d p 2 *1 ^2,1 p 2 -2r I ^2 ? d p 2 

Tl p 2 ^ 2,0 * ri - 1 Tg Ap 2 ^ “ n - 1 Mg Ap 2,o 


1 (bio) 


24 


NA.CA 1H 4021 


Thus, during the time interval between passage of the shock and passage 


Ap. 


of the contact surface, 


2,d 


Au 


remains at about 1, 


2,d 


varies from 1 


Ap 0 Auu, _ 

Ap„ , p„ 2,o AT 2 ^ p 2 2 - > ° 

to 0, 2 — ^ varies from 0 to 1 /y 1 , rp ' ^ varies from 1 to 


n - 1 AM 2 d p 2 

, and — 


varies from 0 to 


r 1 - 1 


In interpreting 


Ti ' % Ap 2?0 — — — 2r 1 

these results, recall that Ap 2 Q and Au 2 ^ Q are negative. It is seen 

that there is only a slight increase in AMg d during the course of an 
experiment . 


Perturbations at Time t for x 



u 2 x 

in Range — — r ^ 1 

u s u s t 

When the variations with 
x at a fixed time t are con 
sidered, subscript 2,o repre 
sents conditions directly be- 
hind the shock at a fixed time 
under consideration. The per- 
turbation equations have the 
same form as in the previous 
section except for As 2 d /c v , 

which is now written 


^2 ,d p 2 
c v Ap 2,o 


(r x - i)(m^ - 1) 2 
M \ [(r x - i)K? + 2] 


l-n 



which follows from 


(AP 2 ,o) . 

x.,t 


A'°A 


Ap 


2,o 


l-n 



(Bll) 


This modification is in agreement with the discussion associated with 
equation (B5), which indicates that a solution for a fixed x can be 
converted to a solution for a fixed t by multiplying the former by 
l-n 
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Figure 1. - Characteristic lines contributing to pressure and velocity per- 
turbations at point d(x,t) behind shock. 
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Figure 2. - Axial variation of pressure perturbation in region 2 at fixed 

time. Air-air shock tube, T-. = T, = T * 520° R. 
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Figure 2. 


Concluded. Axial variation of pressure perturbation in region 2 at fixed time. Air-air shock tube, 

T = 520° R. 
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(a) Laminar boundary layer. 

Figure 3. - Axial variation of velocity perturbation In region 2 at fixed time. Air-air shock tube 
T-. - T, = T = 520° R. 
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Figure 3. - Concluded. Axial variation of velocity perturbation in region 2 at fixed time. Air-air shock 
tube, T x = T 4 = T w = 520° R. 
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(a) Laminar boundary layer. 

Figure 4. - Attenuation in air-air shock tube (ref. 1, fig. 5). T-, = T 4 = 520° R; Prandtl number, 0.70; ratio of specific 

heats, 1.4. 
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(b) Turbulent boundary layer. 

Figure 4. - Concluded. Attenuation in air-air shock tube (ref. 1, fig. 5). = T 4 = 520° R; Prandtl number, 0.70; ratio 

of specific heats, 1.4. 
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(u s t/x) - 1 
(a) Laminar boundary layer. 

Figure 5. - Time variation of pressure perturbation at fixed station x in region 2. Air-air shock tube, 
T x = T 4 = T w = 520° R. 
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(b) Turbulent boundary layer. 


Figure 5. - Concluded. Time variation of pressure perturbation at fixed station x 
in region 2. Air-air shock tube, T^ a T 4 * T w ■ 520° R. 
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Figure 6. - Time variation of velocity perturbation at fixed station x in region 2. 
Air-air shock tube, * T4 =* T w = 520° R. 
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Figure 6. - Concluded. Time variation of velocity perturbation at fixed station x 


in region 2. Air-air shock tube, T 1 


T w = 520° H. 
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longitudinal distance from diaphragm, 8 feet; hydraulic diameter, l/7 foot. 
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(a) Laminar boundary layer. 


Figure 8. - Pressure perturbations (at 
behind weak waves (eq. 12(a)). y^ - 


fixed instant of time) 
T 4 = 1.4; T 1= T 4 . 
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(b) Turbulent boundary layer. 


Figure 8. - Concluded. Pressure perturbations (at fixed 
instant of time) behind weak waves (eq. 12(a)). 

Tq = I 4 = 1*4; Tj = T^ . 
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(a) Laminar boundary layer. 


Figure 9. - Velocity perturbations (at 
behind weak waves (eq. 12(b)). T]_ = 


T, = T 


4 * 


fixed instant of time) 
T 4 = 1.4; 
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(b) Turbulent boundary layer. 

Figure 9 . - Concluded . Velocity perturbations ( at fixed 
instant of time) behind weak waves (eq. 12(b)). 

Tl = r 4 = 1.4; T x = T 4 . 
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Figure 10. - Time variation of pressure perturbations (at fixed station) behind weak waves (eq. 12(a)) 
r ± = r 4 = 1-4; T 1 = t 4 - 


NACA TN 4021 


1.0 


A P2,a 

Ap 2,o 



(b) Turbulent boundary layer. 


Figure 10 - Concluded. Time variation of pressure perturbations (at fixed station) behind weak waves 
(eq. 12(a)). r 1 = T 4 = 1.4; = T 4> 
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(a) Laminar boundary layer. 


Figure 11. - Time variation of velocity perturbations (at fixed station) behind weak waves (eq. 12(b)). 
Ti = T 4 = 1.4; = T 4 . 
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(b) Turbulent boundary layer. 

Figure 11. - Concluded. Time variation of velocity perturbations (at fixed station) behind weak waves 
(eq. 12(b)). r x = r 4 = 1.4) 
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